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Abstract. - Lagrangian chaos is experimentally investigated in a convective flow by means 
of Particle Tracking Velocimetry. The Finite Size Lyapunov Exponent analysis is applied to 
quantify dispersion properties at different scales. In the range of parameters of the experiment, 
Lagrangian motion is found to be chaotic. Moreover, the Lyapunov depends on the Rayleigh 
number as TZa^^^. A simple dimensional argument for explaining the observed power law scaling 
is proposed. 



The investigation of transport and mixing of passive tracers is of fundamental importance for 
many geophysical and engineering applications [Q, |^. It is now well established, and confirmed 
by several numerical j^] and experimental Q evidence, that even in very simple Eulerian flow 
(i.e. laminar flow) the motion of Lagrangian tracers can be very complex due to Lagrangian 
Chaos 1^, In such a situation, diffusion may be of little relevance for transport which is, on 
the contrary, strongly enhanced because of chaotic advection ^. 

In this Letter, we address the problem of quantifying the dispersion of passive tracers in a 
relatively simple convective flow at various Rayleigh numbers TZa. By applying the Finite Size 
Lyapunov Exponent (FSLE) |Q analysis to the Lagrangian trajectories obtained from Particle 
Tracking Velocimetry (PTV) technique, we are able to estimate the dispersion properties at 
different scales. We flnd a clear power law dependence of the Lagrangian Lyapunov exponent 
on the Rayleigh number. This dependence is explained by a dimensional argument which 
excludes a role of diffusion in the dispersion process. 

The experiment is performed in a rectangular tank L = 15.0 cm wide, 10.4 cm deep and H — 
6.0 cm height, filled with water. Upper and lower surfaces are kept at constant temperature 
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while the side waU are adiabatic. The convection is generated by an electrical circular heather 
of radius 0.4 cm placed in the mid-line of the tank at 0.4 cm above the lower surface. The 
heather works at constant heat flux Q which is controlled by a feedback on the power supply. 
By changing the heat input we control the Rayleigh number TZa — {g(3QH^) / {olvk), where g is 
the gravitational acceleration, [3 the thermal expansion coefficient, a the thermal conductivity, 
u the kinematic viscosity and n the thermal diffusivity. In the parameters range explored in 
our experiments the flow consists of two main counter-rotating rolls divided by an ascending 
thermal plume above the heat source The upper end of the plume oscillates horizontally 
almost periodically with a frequency which depends on the Rayleigh number. 

Lagrangian trajectories are obtained by PTV technique ||^. The fluid is seeded with a large 
number of small (50 /im in diameter), non-buoyant particles. The vertical plane in the middle 
of the tank and orthogonal to the heat source is illuminated by a thin laser light sheet. Single 
exposure images are taken by a CCD camera and subsequently digitalized at 8.33 Hz rate with 
a 752 X 576 pixels resolution. Trajectories are then identified as time ordered series of particle 
locations. 

Each run lasts for 2700 s, corresponding to 22500 frames. Typically 900 particles are 
simultaneously tracked for each frame. In the following we analyze the trajectories obtained 
in 6 different runs with Rayleigh number in the range 6.87 10^ < TZa < 2.17 10^. 

The analysis of Lagrangian data have been done with the Finite Size Lyapunov Exponent 
tool which, introduced in the context of the predictability problem turbulence Q , has already 



been demonstrated very efficient for the characterization of dispersion in bounded domains 1 10 
and in the treatment of experimental data Let us recall the basic ideas on the FSLE; more 
details can be found in [Q, |l^ . The idea of the Finite Size Lyapunov Exponent is to generalize 
the Lyapunov exponent, which measures the average rate of divergence of two infinitesimally 
close trajectories, to finite separations. To this aim, we fix a set of thresholds i?„ = i?op" 
(n = 0, . . . , N) and we consider, at each time i, new couples of trajectories Xi(t), X2(i) found 
at separation R{t) = |xi(i) — X2(t)| < Rq. We follow the evolution of the trajectories and 
compute the "doubling time" Tp(Rn) it take for the separation to grow from scale i?„ up to 
Rn+i = pRn- Of course it must be p > 1, but we take p not too large in order to avoid 
contributions from different scales. 

After performing a large number of doubling time experiments (over the possible different 
couples in the run) we average the doubling time at each scale R from which we define the 
Finite Size Lyapunov Exponent 

where <[...] >e indicates the average performed on the doubling time experiments (see Refs. 
for further details). The FSLE is a generalization of the Lyapunov exponent A [0, in 
the sense that 

limA(i?)-A, (2) 

-fl — >u 

physically speaking \{R) ~ A for R < Ie, where Ie is the smallest Eulerian characteristic 
length. For larger values of R, X{R) gives information on the mechanism governing the 
dispersion at scale R. For example, in the case of standard diffusion, on the scales in which 
diffusion establishes, one has 

X{R)~D/R^, (3) 



where D is diffusion coefficient. 
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The use of the FSLE is particularly useful for studying the dispersion properties of passive 
tracers in closed basins [|lO|, and, therefore, fits very well with the considered flow. In such a 
situation, asymptotic regimes like diffusion might never been reached due to the presence 
of boundaries. Denoting by Rmax the average couple separation in the asymptotic uniform 
distribution, it has been found that for a large class of system, for R close to Rmax, fits 
the following universal behavior 

A(fi)^l -^— (4) 

Tfl R 

where tr has the physical meaning of the characteristic time of relaxation to the uniform 
distribution. Equation can be obtained assuming an exponential relaxation of tracers' 
concentration on the uniform distribution [ p^ . 

We applied the FSLE analysis to the Lagrangian trajectories experimentally obtained. In 
order to increases the statistics at large separations i?, we have computed A(i?) for different 
values of the smallest scale i?o (i?o = 0.4 cm, 0.6 cm, 0.8 cm). The threshold ratio is p = 1.2 
for all the analysis. For the results presented below we use H — 6 cm as unit length and the 
diffusive time — H^/k ~ 25000 s as unit time. 

In Figure 1 we show the X{R) versus R computed for the run at TZa = 2.39 10®. The first 
important result is the convergence of A(i?) to the constant value A ~ 3100 at small R. 
This corresponds to an exponential divergence of close trajectories, i.e. a direct evidence of 
Lagrangian chaos. The large value of the Lyapunov exponent (in unit of inverse diffusive time) 
indicates that chaotic advection is the main mechanism for particle dispersion at small scales. 

For larger separation A(i?) drops to smaller values, indicating a slowing down in the 
separation growth. This is quantitatively well described by the saturation regime (^. The 
collapse of the curves at different Rq confirms that sufficient high statistics is reached even at 
these large scales. Fluctuation among different Rq curves can be taken as an estimation of the 
error for A(i?). By fitting the large scale behavior of X{R) with we obtain Rmax — 1.9 7J 
and Tn ~ 8.0 lO"'*^^- Thus also for the late stage of relaxation to the uniform distribution, 
diffusion seems to play a marginal role. The characteristic Eulerian scale in the flow Ie can 
be estimated by the end of the exponential regime (plateau A(i?) ~ A) at which the non linear 
effects start to dominate. We find Ie — 0.5 i?, indeed not too far from the saturation value. 
In this condition there is no room for the development of a diffusive regime [|l0|, as Fig. 1 
clearly shows. In addition, let us mention that by comparing A(i?) computed at different TZa, 
we find that both the characteristic scales Ie and the saturation scale Rmax are independent 
on the Rayleigh number. 

In order to explore the dependence of the Lagrangian statistics on the Eulerian charac- 
teristics, we have performed the FSLE analysis for Rayleigh number varying over more than 
one order of magnitude. The dependence of the Lagrangian Lyapunov exponent A (computed 
from the plateau of A(i?) at small R) on TZa is shown in Figure 2. A clear scaling is observed, 
indicating a power law dependence 

A - TZa^ (5) 

with 7 = 0.51 ± 0.02. 

It is worth noting that because of the geometry of our experiment, the flow shows qual- 
itatively the same pattern for the whole range of TZa explored. This is conflrmed by the 
independence of Ie and Rmax on TZa as discussed above. In these conditions, the scaling of A 
on TZa can be supported by the following dimensional argument. The equations of motion in 
the Boussinesq approximation and made non-dimensional in terms of H and and rescaling 
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the temperature fluctuations T with the typical temperature difference AT, are [1^ 



1 



dUa dUa d 

dt dx/3 dXa^ 



d^Ua 

- TlaTza (6) 



dT dT d^T 

where Pr = u/k is the Prandtl number (which is kept constant in our experiments). It is now 
easy to verify that performing the following rescaling: 

ua Aua , t , TZa ~* K^Ua , (8) 



where A is an arbitrary factor, equations (|6|-|7|) remain unchanged, a part the diffusive terms. 
This means that if one can neglect the diffusive effects (as it is suggested by the previous 
results), the Boussinesq equations are invariant with respect the rescaling (^). Let us stress 
that this rescaling do not involve neither the space (as it is suggested by above discussion) nor 
the Prandtl number (because in our experiments we change Q, keeping both v and k constant). 
The consequence of the Eulerian scaling invariance on the Lagrangian motion, governed by 

^=u(x(t),t), (9) 

is that Lagrangian trajectories are independent on the Rayleigh number. The Lyapunov 
exponent, which is dimensionally the inverse of a time, rescales with TZa^^^, according to 
the result shown in Figure 2. 

The dimensional argument (^||) implies a rescaling of all the FSLE, not only the linear 
part A, i.e. that X{R)/TZa^'^ is a i?a-independent function. In Figure I we plot the FSLE 
compensated with TZa}^'^ for the different runs. The collapse is indeed rather good confirming 
the validity of the dimensional argument and of having neglected the diffusive terms. 
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Fig. 1. - X{R) versus R for different initial 

thresholds Rq = 0.067 (o) , 0.1 i?, (A), 0.13 H (v) at Ua = 2.39 10**. The straight line 
is the Lyapunov exponent A = 3100 ± 200 and the curve is the saturation regime (||) with 
TR = 8.010~4t, and R„,ax = 1-9 i/. 

Fig. 2. - Lagrangian Lyapunov exponent dependence on the Rayleigh number TZa. The errors 
are estimated by the fluctuations at different initial Rq- The line is the best fit A ^ TZa'^'^^. 

Figure captions. - 



Fig. 3. - Data collapse of X{R) at different TZa rescaled with TZa^^^. 
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